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ABSTRACT 

Context. Methanol masers are associated with young high-mass stars and are an important tool for investigating the process of massive 
star formation. 

Aims. The recently discovered methanol maser ring in G23. 657-00. 127 provides an excellent "laboratory" for a detailed study of the 
nature and physical origin of methanol maser emission, as well as parallax and proper motion measurements. 

Methods. Multi-epoch observations of the 12.2 GHz methanol maser line from the ring were conducted using the Very Long Baseline 
Array. Interferometric observations with milliarcsecond resolution enabled us to track single maser spots in great detail over a period 
of 2 years. 

Results. We have determined the trigonometric parallax of G23. 657-00. 127 to be 0.313+0.039 mas, giving a distance of 3.19^3j kpc. 
The proper motion of the source indicates that it is moving with the same circular velocity as the LSR, but it shows a large peculiar 
motion of « 35 km s -1 toward the Galactic center. 
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1. Introduction In this paper we report multi-epoch observations of 
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G23.657-00.127 at 12.2 GHz using the NRAfjQ Very Long 
Two spectral lines of methanol masers, 6.7 GHz and 12.2 GHz, Baseline Array (VLBA). We focus on the measurement of the 
are strongly associated with young high-mass stars or protostars source parallax and proper motion, providing a reliable distance, 
and enable one to investigate the stellar environment on scales Maser variability and internal proper motion studies for both the 



of 100-1000 AU in great detail (Menten [199U Moscadelli et al. 6.7 and 12.2 GHz lines will be presented in a separate paper. 
■ 119991 Minier et al. 2000). Interferometric studies of methanol 
masers have been carried out for more than ten years in order to 

constrain how and where massive stars are being born (Norris et 2. Observations and data analysis 

al.[T998l Phillips et al. [19981 Walsh et al. [19981 Dodson et al. 

120041 1. We started interferometric studies of the 6.7 GHz masers 2. 1 . Observations 

discovered in the Torun unbiased survey towards the Galactic „ , . ^ . . . .. . „. 

, ,„ i.i i ^rvrv-it Our first interlerometric observations or the methanol maser 

plane (Szymczak etal.[2002*. fa G23 65? _ 00 12? at 12178 . 597M Hz using the NRAO Very 

One source, G23.657-00.127, was imaged during the Large Baseline Array (VLBA) were carried out on 2005 

European VLBI Network (EVN) session in 2004 when eight November 28 (epoch 1). After successfuly detecting and imag- 

antennas were working for the first time at 5 cm wavelength, ing the maser we made five additional observations between 

These observations revealed a new class of sp herica lly symmet- 2006 June 10 and 2007 June 19 (see details in Table [TJ. Ten 

ric methanol maser sources (Bartkiewicz et al. 120051 . The "ring" antennas were used during all sessions but only nine anten- 

of masers in G23. 657-00. 127 has a mean radius of 127 mas. Its nas were working in epoch 3, as Hancock was not available 

circular geometry makes it a unique laboratory for detailed study then. Each observation lasted 10 hours and we used one 8 MHz 

of a methanol maser in an isolated massive star forming region, band in dual circular polarization mode. The band was cen- 

However, the morphology and velocity signature of the maser tered at the Local Standard of Rest (LSR) velocity of 83 km s _1 , 

spots at 6.7 GHz did not allow us to determine unambiguously 



the origin of the ring-like structure. Therefore we started VLBI i The National Radio Astronomy Observatory is operated by 
observations of the 12.2 GHz methanol maser line following its Associated Universities, Inc., under a cooperative agreement with the 
detection with the Torun 32 m antenna. National Science Foundation. 
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Table 1. Details of the VLBA observations 



Epoch 



Date 


J 1825-0737 


Beam 


Image rms 




Flux density 


HPBW; P.A. 


per channel 




(mJy) 


(mas 2 ); (°) 


(mJy beam -1 


2005 Nov 28 


166 


3.7x1.2; -13 


10 


2006 Jun 10 


204 


2.3x1.1; -5 


7 


2006 Sep 21 


208 


2.2x1.6; -3 


7 


2006 Dec 14 


219 


2.5x0.9; -12 


7 


2007 Mar 21 


278 


2.2x1.3; +1 


7 


2007 Jun 19 


275 


2.3x1.1; -6 


7 



which corresponded to the brightest feature of 6.7 GHz emission 
(Bartkiewicz et al. I20051 >. 

We u sed the VLBA calibrator J 1825-0737 (Fomalont et 
al. [2003 ) as the phase-reference source. The flux densities 
of the phase calibrator in all observing sessions are summa- 
rized in Table [TJ The pointing position for G23. 657-00. 127 
was chosen to be near the brightest spots detected at 6.7 GHz: 
RA=18 h 34 m 5 1*565, Dec=-08°18'21'/3045 (J2000). The sepa- 
ration between the phase-reference and maser source was ~2°4. 
For epoch 1 we used a 2.5 min switch cycle between the phase 
reference source (slewing plus on-source time = 1 min) and the 
maser source (slewing plus on-source time =1.5 min). The total 
on-source time for the maser was 5 hours. In addition, we ob- 
served the strong source, J 1824+01 19 (Petrov et al. 120051 1. every 
hour for 2 min in order to improve the delay calibration. 3C454.3 
was used as a fringe finder (Ma et al. 1998). 

Starting with the second epoch, we shortened the time de- 
voted to the maser source from 1.5 to 1.0 min, in order to ob- 
tain better phase transfer. We also included a second extragalac- 
tic background source, J1833-0855 (Xu et al.2006), alternating 
with 16 min cycles of rapid switching between Jl 825-0737 and 
J1833-0855 and between J1825-0737 and G23.657-00.127. 
Unfortunately, J 1833-0855 was not useful for astrometry, since 
the source was weak, partially resolved and not detected in all 
epochs. The total on-source time for the maser was 2.3 hours for 
epochs 2-6. 

Starting with epoch 2, we also performed geodetic-like ob- 
servations using 39 ICRF quasars whose positions were known 
to better than 1 mas. Typically 17 quasars were observed over 
a timespan of 40 min, with a frequency setup involving eight 
4 MHz bands at left circular polarization that spanned a fre- 
quency range of 470 MHz. Three of these blocks were placed 
at the beginning, in middle and end of these five sessions. This 
enabled us to estimate tropospheric zenith delay errors and clock 
offsets from the correlator model (see Reid & Brunthaler 120041 
and Brunthaler et al. 2005 for a detailed discussion). 

The data were correlated at the VLBA correlator in Socorro, 
NM. The phase-referencing data were correlated with 1024 
spectral channels, corresponding to a channel spacing of 
0.193 km s _1 . The geodetic block data, containing the observa- 
tions of the ICRF quasars, were correlated with 16 spectral chan- 
nels. 



2.2. Calibration and imaging 

Most of the calibration and data reduction were carried out 
with standard procedures for spectra-line observations using the 
Astronomical Image Processing System (AIPS). First, we ap- 
plied the latest values of the Earth's orientation parameters and 



corrected for effects of the changing antenna parallactic angles. 
.We removed atmospheric zenith delay errors using the geodetic- 
like observations of the quasars. Then we calibrated the visibility 
amplitudes using measured antenna gains and system tempera- 
lures and corrected for voltage offset in the samplers. Next, we 
corrected for electronic delay and phase differences among the 
IF bands using data from J 1824+0 199 and applied the phase- 
referencing technique using J1825-0737. Finally the interfer- 
ometer spectra for G23.657-00.127 were shifted to keep the 
LSR velocity of 83 km s -1 centered in the band. Because, in 
this paper, we are concentrating on the astrometry, we do not 
consider here the results of self-calibration of the maser. A de- 
Tailed disscussion of nature of the 12.2 GHz emission towards 
G23. 657-00. 127 will be presented in a forthcoming paper. 

We imaged an area of 0.4x0.4 arcsec 2 with four maps, each 
covering (0.2x0.2 arcsec 2 ). Each map contained 1024x1024 
pixels of size 0.2 mas and natural weighting was applied during 
mapping. We searched for emission over the Vlsr range from 
72 to 92kms -1 . The resulting beam sizes and rms noise levels 
(lcr) in line-free channels for each session are summarized in 
TableQ] 

The positions of all maser spots (above the level of 5<x on 
the individual channel maps) were determined by fitting ellipti- 
cal Gaussian models (AIPS task JMFIT). The formal errors of 
the fitting were typically 0.01-0.1 mas in right ascension and 
0.02-0. 15 mas in declination, depending on maser strength and 
structure. 



3. Results 

3.1. 12.2 GHz emission 

In epoch 1, methanol maser emission at 12.2 GHz towards 
G23.657-00.127 was detected in a range from 77.41 to 
86.84 km s _1 , which is similar to the range over which 
6.7 GHz methanol maser emission was seen (77.0-87.8 km s , 
Bartkiewicz et al. 2003). In total we measured 34 maser 
spots in the individual channel maps. The strongest spot 
(648 mJy beam -1 ) appeared at Vlsr = 82.60 km s _1 . The dis- 
tribution of 12.2 GHz emission closely follows the ring seen in 
the 6.7 GHz line. Indeed, all 12.2 GHz spots have 6.7 GHz coun- 
terparts but the reverse is not the case. 

During the next five observations emission was seen in the 
same velocity range, but more spots were detected above the 5<x 
limit on the individual channel maps. There were 68, 86, 79, 68 
and 85 spots at epoch 2, epoch 3, epoch 4, epoch 5 and epoch 6, 
respectively. The brightest spot showed a constant LSR velocity 
(82.60kms -1 ) but variable peak brightnesses of 619, 691, 647, 
817 and 911 mJy beam -1 at epochs 2 through 6. As for epoch 1, 
all 12.2 GHz maser spots at all five epochs had 6.7 GHz counter- 
parts. The circularly symmetric morphology of maser emission 
was clearly seen in all five epochs. In Fig. QJ we present total 
intensity map of 12.2 GHz emission from the epoch 6 data. 

3.2. Parallax measurements 

In order to determine the source parallax, we searched for spots 
which were detected in all six epochs, i.e., coinciding in veloc- 
ity and position relative to the other maser spots. We found 19 
such spots and modeled their change in position over the six 
epochs as due to the effects of annual parallax and proper mo- 
tions (see the example in Fig. [2). Since relative position mea- 
surements are usually dominated by systematic errors, owing to 
uncompensated atmospheric delays, formal position errors are 
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Fig. 1. Top: Total intensity ("zeroth" moment) map of 12.2 GHz methanol maser from G23. 657-00. 127 from 2007 June 19 using 
the VLBA. The color scale (in the electronic version) varies linearly from to 180 Jy beam 1 ms _1 . The 19 crosses represent spots 
that were chosen for the parallax estimation (Table|2). Bottom: The enlarged views of the masing regions. The contours correspond 
to the value of 10% of the peak (i.e., 38.6 Jy beam 1 ms _I ). The beam is indicated by the ellipses in the bottom left-hand corners. 



unrealistically small leading to rather large x 2 per degree of free- 
dom values (~ 5) for the fit. In order to account for systematic 
position errors and obtain reasonable formal fitting uncertainties, 
we added "error floors" of 80 and 160 //as in right ascension and 
declination, respectively, in quadrature to the formal fitting er- 
rors. This resulted in x 1 per degree of freedom values of ~ 1 for 
both the right ascension and declination data. 

In Table |2] we list all 19 spots used in the fitting and summa- 
rize the results of the parallax and proper motion fitting individu- 
ally for each sequence of spots: the mean parallax is 0.3 12 mas, 
with a standard error of the mean of O.OlOmas (standard de- 
viation of 0.042 mas). A combined fit for all spots with a sin- 



gle parallax, but individual proper motions and position offsets 
for each spot (since all the maser spots are at the same dis- 
tance within the measurement accuracy), gives a similar result 
of 7r=0.313±0.015mas. 

In Fig.|3]the combined fit for the parallax with data from four 
representative spots (spots 1-4 from Table |2]l is presented. The 
positions of the masers after removing the proper motion and 
position offset are shown for clarity. Combining the results of 
several maser spots can lead to underestimated parallax uncer- 
tainty, since the position measurements relative to a background 
source may not be independent. Random errors (e.g., from the 
position fits due to the finite signal to noise ratios and possible 
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Fig. 2. Left: Six epoch positions of spot 19 (from Table |2} on the plane of the sky with the best parallax and proper motion fit 
superposed. Open circles indicate the location of the model value at the dates of observation. Right: Position of spot 19 versus time 
with the best parallax and proper motion fit in right ascension (filled circles, solid line) and declination (squares, dashed line). 
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Fig. 3. Position offsets of four spots (filled circles) in right ascension (top) and declination (bottom) after removing the proper 
motion, superposed on the combined parallax fit model for all maser spots in Table[2] The data points are slighlty shifted in time for 
clarity. Also shown are the average positions of all spots (black triangles). 



maser spot structure) are probably not correlated among different 
maser spots. However, systematic errors (e.g., caused by residual 
atmospheric delay errors) will affect all maser spots in one epoch 
in a very similar way. Thus, the formal error in the combined fit 
will underestimate the true error. 

The most conservative approach in the uncertainty estima- 
tion would be to assume that the errors are 100% correlated 



and multiply the formal error with V^V, where N is the num- 
ber of maser spots. This would give a parallax uncertainty of 
±0.015 a/19 = 0.065 mas. To estimate the effect of the system- 
atic errors on our parallax measurement, we have calculated the 
average position of all maser spots in each epoch after remov- 
ing their fitted proper motions and position offsets. The aver- 
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Table 2. Spots which appeared at all six epochs. Presented coordinates (columns three and four) are relative to the phase center 
(RA=18 h 34 m 5 1^565, Dec=-08°18'21'.'3045) at epoch 2. Column five lists the parallax estimates. Columns six and seven give the 
motion on the plane of the sky along the right ascension and declination, respectively, ura is the true RA motion in masyr 1 and 
includes the 15 cos(Dec) factor. 



Spot 


Vlsr 


ARA 


ADec 


Parallax 


Proper motion 










n 




/^Dec 




(km s- 1 ) 


(mas) 


(mas) 


(mas) 


(mas yr ) 


(nicis yr ) 


1 


86.837 


-38.090 


16.134 


0.299±0.071 


— 1 .ZU±U. 1U 


7 00-1-0 1 f\ 
— Z.VU+U. 10 


2 


86.452 


-36.343 


13.110 


0.330±0.100 


— l.Zo+U.l J 


7 C7-I-0 7 1 
— Z.o / +U.Z1 


3 


86.260 


-36.149 


12.916 


0.306±0.095 


— 1 .Zu+U. 1Z 


7 QQ-i-O 1 f\ 

— z. yy±u. id 


4 


84.335 


-152.281 


2.371 


0.332±0.070 


1 zL^-t-O OQ 


1. 1 7_i_0 1 7 
— J. 1Z+U. 1Z 


5 


84.143 


-152.105 


2.719 


0.273±0.057 


i t"7_. o m 

— L.J /+U.U / 


"3 07_i_0 1 o 

— o.u / ±u. iu 


6 


84.143 


-71.635 


30.214 


0.269±0.058 


1 9^-t-O 07 
— l.ZO+U.U / 


7 QQ_i_0 1 7 
— Z. J J±V. 1Z 


7 


83.950 


-17.850 


-20.833 


0.336±0.054 


_i ^+0 Ofi 

1 . .J .J LEV/. \J\J 


07+0 OQ 


8 


83.950 


-71.029 


30.442 


0.245±0.052 


-1.43+0.07 


-3.02+0.11 


9 


83.373 


-65.264 


31.895 


0.373±0.070 


i i o_i_n i o 
— i . i y±u. iu 


-1 OO-I-O 1 (\ 

— j.UU+U. 10 


10 


83.180 


-65.064 


32.004 


0.311±0.070 


-1.32+0.08 


-3.03+0.13 


11 


82.988 


-19.591 


-5.390 


0.290±0.064 


-1.37+0.08 


-2.81+0.11 


12 


82.603 


-14.350 


-14.133 


0.277±0.099 


-1.33+0.13 


-3.15+0.17 


13 


82.603 


-18.516 


-1.136 


0.363±0.085 


-1.30+0.11 


-2.75+0.15 


14 


81.833 


-20.069 


18.417 


0.301±0.056 


-1.27+0.07 


-2.87+0.11 


15 


80.486 


15.676 


-171.016 


0.406±0.096 


-1.43+0.11 


-2.87+0.19 


16 


80.294 


16.118 


-171.487 


0.355±0.099 


-1.35+0.13 


-2.96+0.18 


17 


80.101 


16.429 


-171.892 


0.311±0.074 


-1.35+0.10 


-2.91+0.14 


18 


79.909 


16.871 


-172.332 


0.249±0.065 


-1.26+0.08 


-2.85+0.12 


19 


79.716 


17.326 


-172.979 


0.304±0.053 


-1.38+0.07 


-2.92+0.11 


Averas 


'ing fits 






0.312±0.010 


-1.32+0.02 


-2.96+0.03 



Combined fit 0.313+0.015 
Averaging data 0.313+0.039 



age positions are shown as open triangles in Fig. [3] Such av- 
eraging should reduce the random error, but leave the system- 
atic error unchanged. A parallax fit to the averaged data points 
yielded a value of 0.313±0.039mas, which we adopt for the 
parallax of G23.657-00.127. This corresponds to a distance to 
G23.657-00. 127 of 3.19+°^ kpc. 

The average proper motion of all spots and the standard er- 
ror of the mean is -1.32+0.02 masyr -1 in right ascension and 
-2.96±0.03 masyr -1 in declination. For the measured distance 
of 3.19 kpc, these values correspond to: -20.0 ± 0.3 km s _I 
and -44.7 ± 0.5 km s _1 in right ascension and declination, re- 
spectively. All spots have very similar proper motions (within 
0.2 masyr 1 of the average) (Table|2]l. Hence, internal motions 
are small, as is commonly found for 12.2 GHz methanol masers 
(Moscadelli et al. 2002) and, thus, the average proper motion 
closely represents the motion of the young star that excites the 
masers. 

4. Discussion 

Using our measured parallax and proper motion of the 
source G23. 657-00. 127, we can estimate its 3-dimensional 
motion in the Galaxy. Assuming IAU recommended values 
for the distance of the Sun from the Galactic center of 
Ro=8.5±0.5kpc and the circular rotation speed of the LSR of 
#o=220±10 kms -1 , the Hipparcos Solar Motion values (Dehnen 
& Binney [19981 1, our distance measurement of 3. 19+0.40 kpc 
and proper motion of -1.32±0.02mas yr~' in right ascension 
and -2.96+0.03 mas yr 1 in declination, an LSR velocity of 
83+2 km s~', we obtain a Galactocentric motion of: 

U' = 42.3 ±4.0 km s" 1 , 



V = 222.4+ 11. 3 km s _1 , 
W = 4.1 + 1.4kms _1 . 

Here, U' denotes the radial velocity component toward the 
Galactic Center, V is the velocity in direction of Galactic ro- 
tation, and W is the velocity toward the North Galactic Pole. 
Hence, G23. 657-00. 127 is moving with nearly the same rota- 
tion speed as the LSR and slowly out of the plane of the Galaxy. 
However, it has a very large peculiar motion towards the Galactic 
center. This large peculiar motion explains why the source is 
closer than the near kinematic distance of » 5 kpc. 

If we use a different rotation model of the Milky Way with 
Ro=8.0±0.5 kpc and #o=236±l kms -1 which is consistent with 
the measured proper motion of the Galactic center Sgr A* (Reid 
& Brunthaler r20041 >. we get: 

U' = 33.6 ± 3.410ns" 1 , 

V = 238.8 ±5.8 km s" 1 , 
W = 4.1 ± l^kms -1 . 

Here, we find again a similar rotation speed as the LSR and a 
large motion toward the Galactic center, demonstrating that this 
result is not sensitive to the Galactic rotation model. 

With a Galactocentric distance of 5.7 kpc (assuming 
R =8.5±0.5kpc), G23.657- 00.127 is located close to the 
Galactic bar (Benjamin et al. |20051 l. Hence, the large peculiar 
motion might be caused by the gravitational potential of the cen- 
tral bar. Alternatively, the large peculiar motion could be caused 
by an interaction between a density wave of a spiral arm with 
the molecular could. More sources in this region are needed to 
investigate this possibility in more detail. 
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With our new measurement of the distance, we are able to 
verify our earlier calculations of the linear size of the 6.7 GHz 
ring as well as the bolometric luminosity using the mid- and far- 
infrared emission (Bartkiewicz et al. 120051 1. Assuming the mean 
radius of 127 mas the linear size is 405 AU. Applying the for- 
mula of Walsh et al. ( 119971 their Eq. (3)), we get a bolometric 
luminosity of < 11.4 x 10 3 Lq, which could be provided by a 
single B0.5 ZAMS type star (Panagia [1973"j l. 

5. Conclusions 

Using VLBA observations of the 12.2 GHz methanol masers in 
G23. 657-00. 127, we have determined an accurate distance of 
3. 19^Q35 kpc to the maser ring. At this distance, the radius of the 
ring corresponds to a linear size of 405 AU. We have also mea- 
sured the proper motion of the source, which we use to deter- 
mine its 3-dimensional motion in the Galaxy. The source moves 
with a similar rotation speed as the LSR, but shows a high pe- 
culiar motion of 30-40 kms -1 toward the Galactic center. This 
large peculiar motion might be induced by interactions with the 
Galactic bar and is responsible for the large error in the kine- 
matic distance. 
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